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ار ۱۱۱۱۹ 


عصتا0ع211 وعووععم۳م «صفدظ 60امصتصصناه بم 0عنام‌حزه فقط صفاماميم لمممتعممصصتلعمصن مه قه قمع مطا عصا سم 
,60۵۲۵11۷ ۲۱۷۵۲۵۰ طا راتلدتان ما۷۵ عصتاهع]ه 5عمام؟ عطا ۵۶ مج مه 0۳۵۵۵8865 500۲۵26 1۱۷۵۲۰ عط صا تعتفصه] اتصنلمه 
-عع۲ط] مصح لفصمتعصمصتل-۱۳۷۵ 0۶ مویا مطا معصتلومی ع20 تقصمتاه وهای گم صمتاح‌اتصصنا عمط رتماعع؟ لمتاصهاوطاتاه ۵ وه 
ماو ۱۷۲۵2۲ .0619مصظ اقومتفممصصتلعصی 0ع7ماموصصه ب«امل1 9۵۲۵ 0۶ عوا مطا م6 عصتع2ع1 رواممصه لهصمتوصمصتل 
معط تعلنقومع مه مل مصح نله ممتاموظ ممتفهم8ز۲-صمتا۷6۵ 0 فص مق 0عقعه 087610060 میم 618مقظ ععتتام‌طمم 
فص ع0نفجم ما ۲۵۵0960م ود نها ممتاموظ مملومموز- ۸0۷6۵۵۱۵۵ ۱۷۲۵01160 رعوممتتام فنطا ۳۲۵۶ .عصمی2 6ع2عماه 
جمعه فقط ممتبامد 2 ,لته وتطا ما اصمته‌گاوم ممتمممعنل هه ما۷۵ ع۷تا0عه عصتصتلمل اه مه اصقصیماه 
6 .۵0619 لقصمتوصعصصنل-عصن صا ممتاهانصته تلمتان تمه ما عصمی2 عصقصیهاگ فطا ۵۶ )0ع11ه عمط امد ما ۳۳0۵۵0960 
0 36 ۵۶ امممصه۷ م1 مطا ۵۶۵۷۵ ما بام200ظ فنطا ص .عصمتاتلجمع لهعنهه معط م6 تعوماه وا مفمت)ها اه ۲۱۷۵۲ 
عصا 0۶ مه لمطمتاهعو-ووم 2۷67226 فطا ۵۶ حمتامد1 معط وه 8860ع۲مه طمتام1۳2 ععم2 اصقصعهاه 2۷6۲2826 مطا رل‌مطامظ 
معط مه هآه۷۵۱۵ ۵11600۷۵ عطا ۶ه ممتاتصتلمل فطل طمنامطا اع0مص تقصمتوصمصصتلعمصی ه صا ۵0«مامصی فص () ۲۱۷۵۲ 
۷۲ ۵۷ 007610۵60 واع0ممط لهصمتعصمحصتل-مصم مطا ۵۶ ۷۵ تاماه8۵ع۲6۵۲ 106۲60قمع ور [0معظ فنط 1 ,اصملهمه طملوومولل 
۸۳۷220 فط) ‏ 10۲ 0ماهعتاو۷6 1 مه 0مطامصظر . 0عومممعم فط 4ص مصمتاقاه: صمتقفهم‌وزنا ممتاهعت لد فطا 
۰ ۱) ۷۵11024 مه میهعطاتلی م) ۱۷۹۵۵ ۷۵۲۵ ۲۱۷۵۲ عط عممله ما02 لمعمتاه وان 


۵۹ 2۱0 واهز۱۷]۸)۵ 

۶ 02960 00۷610060 ۲«الدعصعع هته فصصهتعهعج عناوم آمصمتوصمصصتلمصم املع هه ماطملته27 مطا من 
و۱0۵0 وه طمناهو )0۲مقصعا راتصللده عصتاهعلله دتمامه؟ عمجم ]مه مل رفظ رلع0مصه مملو۵مو۸0۷6۵)08-1 نا مطا 
0 /0 276 ۵۲006۵8565 0696 ,۳/2 2 طمناو ص .عقعته عمهعماه وه م«مص! فعتاتمصمعمع)عظ اعم1مطامعمحه م1۱ هه 
۶ 26010207 6 عووععص1 ما عع0تن صا رطم‌تهعوع فاص مصمتفنوه لهممتومهممصتل عون 2 وه صملمامعج مطا عصتاصهوميم ره 
0 .0۳0008960 18 متام ۵۲60۷۵ ه واممهص لفصمتفمممتل موم ماماهاتهه مه یمام ی صمتافاتصته معط 
۲ 21۷2۳0 127۷ صا0حاوصه صه 022 مطا رافتت فطا 21 ,نله وعهه عطع مج 0مطاممدر 014060 مه 0عومم۵۲۵ عظ عتهمعمم 
همم لمع موم ۲۱۷۵۲5 م1 ,1۷۵ 0۵0-]۷۵۱ 2 وا طمتطا رهم ناه 2 ف 2006160 15 
من مهد اهتدم ففعصطونا0؟ و مصتصصه]۱۷ .2014 ص1 0242 ماطماتهه معط م6 عصتل؟0عع2 0عصلع0 م۵ فعمتاتصمم 
عص 200۱۷ م) تع0ت0 ص .عمتلناای فناماعم ما عصتل20001 ۳۸۵/5 285 ۵8 تلهم مملوتهم‌ونل معط 24 ,0.033 1۵ 0.017 
18 ۲۱۷۶۲ 06 10۶ ۱۱ ۵۶ مبطلج۷ مطا مصتصصماعل ما لقتاصوووع و1 16 روصمتاهنو 018160مصه مطا ا قهم‌به اصقصعهاه 0۶ اهعگاه 
دا زرمع 1ج پرجمامی صرح برد و٩‏ م2 و(ه) کصمجممون ممتوته‌م‌ولل وره) ۲مامع طم1وهم‌فنل ۵۴ کتعامصصح‌ته معتط) ومونا فتاه 
۰ ۱۷۵۵6 ۲۲۷۵۲ ۸۳۷2۵0 بقعت2 ناد مط مد ممتهاه صقطمگ صا فصمتاتاع6۵: 100 صح فاصه اهتانطا ۵۶ ممتاتصتلعل مطا 
فطل ط ر۳۵7۵]01۲6] ,094 ع< و9 ,273.4 < 2 ,0.168 < ۲ ۲۵۹۵66۵۷۵۲ و27 وطعامعصصعتجم 0عاعصصتاوی فطل ۵۶ فمتاله۷* مط 1 
0 3 مه )مهم ممنوممون 00۷۵عه مه ۵۵0و فطا مصتودا فص م۵ فعمتاهعترهن راعل‌مصه 2۵011160 
۷5 ۱ مه عصنل‌مومع0 اصمت‌لهم ممتوهموتل ماصامته۷ عمنم0زموم 24 0مطامصه نله مملوتهم‌وزنا - م۷6۵0 

)۱۳۷۷6۵۰ جع ۳۲2) وصمتاهاه 0 تعطان عطا صد عم ت۷2 ۷۵۲۵ دععصفطل عوقط 1 ,1ع0مجط لقصمتعصع‌صصتلعصم معط ما ۵660و 


رهظ ه تمد راصمجصم‌م‌مصه]۱۷ مصه مطنهءصتمص عماج ۷۷ 1 ۲۳۵16990۴ ماهههووظ 240 1۳۵۵۲۵۵ وتماو۱۷]2 -3 4صه 2 ,1 
۰ 11210 رهطم 1 ,۲۷۵۲۹16 وعع۷]۵02 120121 
(1268.26.1۳ محط ۵) 1تعطم ۱13222 :اتممصظ ممطه مصت‌طمصوعتزمن) -۶) 


۲ نشریه آب و خاک. جلد ۳۶. شماره ۶ بهمن - اسفند ۱۴۰۱ 


تمه وز(۱ 20 وغاباععع1 

فص 240 1150 0 5۵۲165 مصصتا مطا ۵۶ علههم مطا 0عقننهه ۱۱ تمامصصه‌تهم معط معصتعهعتمصا رقالتاوی: ولبتای عتط ی 60قعظ 
6 وموهع۵ع1 طعتامطاله وع1۳ مط صا فص مها م۳2 عط گم مصتصماتمطه م1 .عفهعععع0 ما فصصتا ۲۵۷۵ ۲۱۷۵۲5 
۶ 56۲165 عصصت معط هن تمامصصتحته فنطا 01 611601 مطا رعصم2 اصقصوفاه عطا ۵۶ ممعممتالصا مطا ما مب )صم‌تگومه ط10وومولل 
و عمصصنا عصلفت مه عصتلاه ۵۶ مومله معط باعل لهصمتاهعهوواه مطا معلنا ,۱۵006۵0 ود ممتاقاجهعومع 12660نصصلو فط 
کمم‌تلههه ممنوه1150 611600۷6 مظ) هصلوتا ۵۶ وععمه ۲۳۵ عطع مد [06محظ آهممتفممصصیلمصن فص مصتتهم‌صصمی بزرظ .1۳20162860 
123 - ومتاماه صقطاهگ 26 0مصتصتماعل قه7 ممتاقاصصته مطا صا پرمح‌تتاهعج ها مقهع:1۳6 عط با ماه مضه بجاه۷۵۱۵ 20 
.امه وم لوعتافتاهاه معتطا طه فاصمه‌تگلهمع فطع عصمتتامه ماه - ماله ممتامطم عبت مطا ۲0 
۹ 386 1۳0 حصع1 150 هه 180 ۵۶ وممصم)عن ده ععبزر مط) عصمله فصمتافاه تعطاه ۳۵ 2 فععصقطل 96عظ) رصم)2001 
۵۲ 19 ۵۶ فطاممصه 12 صا ممتاهاجهعصمي 0مساوحعصهر 0صع 0عاهاناصصه مطا رععصقافضا ۲۵۲ .نامع فتطا دمن ما80 
۶150 15)2206 ۵ 26 0.97 0 0.86 ۵۶ () اصمه‌تگههه ممتاهآهترهه معط دمطمهع( فاصمتم‌تگاههه 2۵0تصصنامن مطا راومه ره 
ص ممص6ع0 مصتصنمنصی: مظ1 ا۵۵ 1:27 1۲0۷۵۵ (تره۴) مه متمناوو فصوعع )۲۵0۵ مظ) هه راون مفطا م0 م1 
نموه ۵۶ مه مطا موه ۵ قعامصصهتهن مان 0 امه مطا رها 0عفننو ۵۵ «قمه ممتامصصتاوه ممتاقتاموعجمع فص 

۰ 136 21002 و8صق1 مص فاصنا 


جمتعصل‌صه) 

6 ۵۶ موه همع 2۳۳۷۵۷۵ 2۷۵ مصتهم‌صلومه ۷۵۲ ما ممتاقتاموعوم هم ممتاح‌امصته مه ممتاقتاوع ملهتتاهعظ 
لهم0قصمصصتلعصم مصتفت ممتتهلیصه ‏ تمتان ‏ ما۳2 1۳010۷6 م6 2960 )موم فلط1.. .وعمممللفص . اهامممصممت وم 
هسامح مطا ما تعوماه مصرمععها 0ج عصتام0مهص عط ۵۶ برمهتتاهعه مطا مفهعتمصا مه تعقتن ظا فتم۷ ۱۱ 0مرت۵1[۱/ م1 2۵06 
5 .60]٩6عوبناه‏ م۵ع0 2۷۵ 1۱۷۵۲ عطا عصمله وعصمم2 0620 مه عصتتملزمممع مه هماو ومامو ممتامعتزهه رفصمتاتل‌جمع 
معط 20 رفحعته اصفصعماه مه مممتاهمو وم و۳۷۵۵ 20ظ عطا ۵۶ ممهلتته مطا ۵۶ 169۵ ,2۷6۲220 08 متقطا ٩00۷760‏ 
136 ۵۶ وممتم0صهمصه 20 ومیل ۵۵0 م1010 قهعتع وعفمط]. ,0مموو و1۱۷۵ مه هم فلصموع0 )صمت‌تاهمه مم1وتهمولل 
ح که 06۵0و ملناوعع م1 .مزر مط ۵ ومهته مطا مه فامصصفط ممتاهعتا 02 معط ود متاصمتاه عامهطاه اقطا کمتوم مط 1 
۰ 0.84 220 ,1.27 ,1.78 10 060162865 عم وعقتانو فصقفطظ )۲۵0 فص رفصمتاهاه ۳۵۷ 200 رحامع/1۳ مصقده5 
2۵290 معط ما تمومام (ماتصتلهه رای متطا م6 ممتح‌اجهعصمن فطل 0متقصصتاوی 061مهصه م1 20011160 مصا م7610[ 
۷ 1 لاه ,۵۷۱0۵0 و1 عصتتم0صهممه ۳۱۷۵۲ وه طمناو وقصم2 0620 0۴ 611601 فطا رفصمتاهاه صقطهه مه طامع۱۳۱۷ 1 .0242 
۲ 36 ۲۵ «اتت(0۳0 وا ما فیطل 09۵98۷۵0 ۷28 وعص2۵ عصقصعهاو 0۶ )مهمرصظا مط صد کصمصصه 1۳0۲0۷ اصهعاگتصعلد مظ بحمتاماه 
مصملاه لته ادن مه [2طم1فصمصصتل عمجم صا پرمه‌تتاهعه ماوت 10۲ عفن ۵ وی طمت2عوع7 فقطا ۵۶ فااناوعم م1 ,طتامظظ 

۰ 1 مطمت) نومه احساقط مط مه تعوماه قرعم مط عصلعص9 200 


6 2201۳20/6 ,1۱۷۵ ۸۳۷۵0 رحمتاهمتصتامه تحمامی کح رصمتاهناوه 1۲16۵۵۲5108 -ص۵۵)0 ۸۷ 1۵۲۱۷۵۲۵۵۰ 
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بجلق ۶ شباره من تفت ۲۰۱ این ۳۷۵ ۳۶ 


بررسی اثر رفتار نواحی ماندابی در مدل‌های یک بعدی هیدرودینامیک و انتقال مواد محلول در 
رودخانه‌ها 


امیرحسین منتظری "- سجاد خدام باشی امامی ۲- مهدی مظاهری " 
تاریخ دریافت: ۱۴۰۱/۱۴/۲۰ 


تاریخ پذیرش: ۱۴۰۱/۰۹/۱۶ 


چکیده 


بیان رودخانه به عنوان مسئلةٌ یک بعدی باعث شده بسیاری از عوامل دخیل در انتقال آلاینده در رودخانه ساده شوند و یا در نظر گرفته نشوند. در 
رودخانه‌ها نواحی نگهداشت عاملی اثرگذار در کیفیت آب می‌باشند. موم محدودیت داده‌های مشاهداتی که یکی از عوامل محدود کننده استفاده از 
مدل‌های دوبعدی و سه‌بعدی است باعث استفاده گسترده‌تر از مدل‌های یک‌بعدی شده است. اکثر مدل‌های تجاری موجود بر اساس معادله انتقال- 
انتشار کلاسیک توسعه یافته‌اند و نواحی نگهداشت را در نظر نمی‌گیرند. به این منظور برای مدل‌سازی تأثیرات نواحی ماندابی از سرعت و ضریب 
پراکندگی موّثر به عنوان روشی برای بهبود عملکرد معادلة انتقال- انتشار پیشنهاد شده است. سرعت و ضریب پراکندگی موثر از طریق اصلاح سرعت 
متوسط جریان و ضریب پراکندگی طولی, با نسبت مساحت سطح مقطع ناحية ماندابی به سطح مقطع عرضی رودخانه تعریف می‌شوند. با جای‌گذاری اين 
پارامترها در حل معادلات سنت- ونانت و انتقال- انتشار در مدل‌های یک‌بعدی می‌توان نتایج را بهبود بخشید. در این مطالعه برای اثبات بهبود عملکرد 
مدل از 1۷1116۳11 به عنوان نماینده‌ای از مدل‌های یک بعدی اعمال و رودخانة اروند با آن مورد بررسی قرار گرفت. برای واسنجی و صحت‌سنجی مدل 
از داده‌های مشاهداتی غلظت در طول رودخانه استفاده شد. پارامترهای فاکتور پراکندگی و توان پراکندگی همراه با نسبت مساحت سطح مقطع ناحية 
ماندابی به سطح مقطع عرضی رودخانه برای ایستگاه 56 در ۱۲۳ کیلومتری از القرنه» با الگوریتم بهینه سازی کلونی مورچگان برآورد شد و برای 
ایستگاه‌های ۳۵۷۷7 و 2۷766 با داده‌های غلظت مورد صحت‌سنجی قرار گرفت. در ایستگاه 56120 12۷۵۵0 و ۳2۸۷ سه پارامتر خطای ضریب 
همبستگی» نش- ساتکلیف و جذر میانگین مربعات در مدل اصلاح شده کاهش یافته و به عنوان مثال پارامتر آماری خطای جذر میانگین مربعات به 


ترتیب به ۰۱/۷۸ ۱/۲۷ و ۰/۸۴ کاهش می‌یابد؛ بنابراین اصلاحات حاکی از بهبود نتایج پیش‌بینی طبق مدل مفهومی جدید است. 


واژه‌های کلیدی: اروندرود الگوریتم مورچگان معادلة انتقال- انتشار ناحیة نگهداشت 


مقدمه 

رودخانه‌ها و نهرهای طبیعی منبع اصلی آب آشامیدنی» کشاورزی 
و صنعتی برای سال‌های متمادی بوده‌اند. در سال‌های گذشته حفظ و9 
بهبود کیفیت آب آن‌ها برای محققان و9 دولت‌ها اهمیت بیشتری پی دا 
کرده است. برنامه‌ریزی و مدیریت برای برداشت و استفاده از رودخانه 


۸ ۲ و ۲- به‌ترتیب دانش‌آموخته‌گان کارشناسی ارشد و دانشیار گروه مهندسی و 
مدیریت آب دانشکده کشاورزی دانشگاه تربیت مدرس, تهران ایران 

1۲. ۱00168 ۵) ترعطم ۱۱۵2 ظ :تم 
1 :]10 


(#- نویسنده مسئول: 


مستلزم درک کامل از کیفیت آب و فرآیندهای موثر بر سرنوشت مواد 
آلاینده در رودخانه است. در رودخانه‌های ساحلی یا جزرومدی که 
منتهی به درا یا محیطهای آبی شور هستنده مصب رودخانه محل 
ورود آب‌شور به آب شیرین رودخانه است (2019 ,۵1 ۶ طنحهعدکل. 
در این رودخانه‌هاء جزرومد باعث نفوذ مکرر جبهةٌ آب شور به رودخانه 
شده و عامل اصلی کاهش کیفیت و افزایش مواد محلول در آب 
رودخانه می‌باشد. این اتفاق موجب کاهش کیفیت خاک اطراف 
رودخانه. کاهش کیفیت آب‌های زیرزمینی و پرهزینه شدن برداشت و 
تصفية آب می‌شود (2019 .۵ 6۶ ۳8 بررسی شوری در 
رودخانه‌های جزرومدی به خصوص در نزدیکی مصب به مراتب 
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پیجیده‌تر از شرایط دیگر است به این علت که پدیده‌ی جزومد دائماً 
باعث تغییر حرکت جریان آب شور دریا می‌شود که در هنگام مد به 
سمت بالادست رودخانه و در هنگام جزر به سمت دریااست. 
شاخص‌های گوناگونی موثر بر پیشروی» اختلاط و نفوذ آب‌شور به 
رودخانه هستند که از آن‌ها می‌توان به ارتفاع جزرومد» دبی‌های 
ورودی به رودخانه» توپوگرافی و ناهمگونی‌های مورفولوژی رودخانه 
که به عنوان نواحی نگهداشت شناخته می‌شوند» شیب بستر رودخانه» 
بت باه ی احتطکا کتسم اشارف کرد که تیم اف را شر 
سرنوشت جبهة آب شور می‌گذارند ,طعصن۹ :2005 ,عزنمه5۵۷) 
(2008. این فرایندها غالباً در مصب رودخانه‌های اختلاط کامل با 
طرح مسئله به عنوان یک معادلةٌ یک بعدی ساده می‌شوند 2۶ ۷۷) 
(2022 ,.۸1 0۶ ۹۵2021 :2012 ,./۸. نفوذ جبهه آب شور در رودخانه 
های جزر و مدی توسط دو فرایند انتقال" و انتشار قابل تشریح است. 
در فرایند انتقال جریان آب همراه خود شوری را منتقل می‌کند و 
فرایند پراکندگی" به علت تغییر در پروفیل سرعت اتفاق می‌افتد و 
باعث باز شدن ابر جبهة شوری در رودخانه می‌شود ,2۳72)) 
(2008. 

در سراسر دنیا مطالعات گوناگونی در رابطه با نفوذ آب شور به 
رودخانه‌های منتهی به گستره‌های آبی شور انجام شده است 
و.]6 6 ۳۵۲6۵۲۵ :2020 .]۵ 6۶ ]۳۱۵000 :2019 ,.۵1 6 تحصحاد۳۲) 
(2022. رودخانه مرزی اروند نیز مستننا نبوده و در سال‌های گذشته 
کیفیت آن مورد توجه واقع شده است :2017 .۵1 ۶۶ طداانتده) 
(2022 ممالضا :2019 بامتاه۱۳ظ :2021 ,.۵1 2 ۸1۱-۵652۷01. به 
عنوان مثال گنگ و همکاران (2022 ,.61 61 0082)) به منظور 
تدقیق اثرات تغییرات دبی دامنه جزرومد و باد را بر نفوذ شوری» 
مصب رودخانه ای در کشور چین را بررسی کردند و نتیجه گرفتند که 
اثرات دامنه جزرومد بر میزان نفود شوری به رودخانه موثرتر از دیگر 
فاکتورها است. به علاوه در مطالعاتی به اثرات اهمیت مورفولوژی و 
بستر رودخانه بر کیفیت آب رودخانه‌های جزرومدی که در محرض 
نفوذ آب شور هستند توجه شده است ,0۷0 :2002 ,.۵1 61 2762ظ) 
(2006 .۵1 21 ص۷۷ :2019 ,.آه ۵71 :]۱۷ :2022 

رشد روزافزون نگرانی‌های زیست‌محیطی سبب شده مطالعات 
گسترده‌ای بر روی فرایندهای موثر بر انتقال مواد محلول در 
رودخانه‌ها انجام گیرد 4حه صدرنطع22 :2021 ,.آه ۶1 ظ۲«۷۵) 
(2022 ,:«161206. فرآیندهای نگهداشت یک مفهوم برای تشریح 
اثر مناطقی درون يا نزدیک رودخانه بوده که سرعت آب در آن‌ها 
نسبت به جریان اصلی کمتر است و می‌تواند شامل استخرهای کنار 


1- 8 
2- 0 


رودخانه, اددی‌هاء نواحی هایپورهیک" و نواحی ماندابی" باشد ( ن) 
0 .1 201667 .]۵ 6 2619106112 :2019 ,.۵1 61 ۲۵202 8 
در استفاده از معادلات یک بعدی در پیش‌بینی کیفیت رودخانه‌هاء 
فرایندهای نگهداشت عموماً ساده شده یا در نظر گرفته نمی‌شوند 
(2013 ,.۵1 ۵ تعطعلایک :2015 رعصهن ده 108۵ ). مقطع 
عرضی رودخانه (۸) به دو بخش شامل ناحية ماندابی (۸۰) و ناحية 
جریان (۸۸) تقسیم می‌شود. مناطق ماندابی در نیرها و رودخانه‌ها 
معمولا در پشت تلماسه‌های بستر» در نزدیکی دیواره‌های مقعر در 
پیها روک وبا کنیل در وی رود عانه کین سی تیوه 
(شکل ۱). پهنه‌های ماندابی قسمتی از رودخانه يا نهر هستند که از 
جریان اصلی جدا شده‌اند. اما به طور متوسط در محاسبه مساحت 
مقطع عرضی کل رودخانه به حساب می‌آینده از ین رو مناطق 
ماندابی» سطح مقطع موثر جریان را کاهش می‌دهند ۶ 6۳6)) 
(2022 ,۱۷۲۵22۵۲1 0مح تطمد‌طاصع100 :2021 ,./ه. این اتفاق 
باعث شده سرعت موثر بالاتر از سرعت میانگین اندازه‌گیری شده و یا 
محاسبه شده از طریق روابط ریاضی باشد ,ع08ع6) 24 0ع۹) 
(2001. سرعت موثر به صورت زير محاسبه می‌شود: 


)۱( ۸ - < , ا 


۸7 سطح مقطع عرضی ناحية جریان* | » سرعت نقطه‌ای در 
ناحية جریان 11 آن]]» ۷7 سرعت موّثر جربان در رودخانه [""آ1] و ۸ 
سطح مقطع کل رودخانه [*1] است. در مناطق ماندابی به علت اینکه 
انتقال (جابه‌جایی) پديدة غالب نبوده. سرعت در آن‌ها از سرعت 
متوسط در رودخانه کمتر است و به صفر میل می‌کند. با توجه به این 
۸ > ۸۶ می‌باشد در نتیجه 7 < 72 است. در اینجا [۷ سرعت متوسط 
آب در رودخانه [1/171] تعریف می‌شود که از تقسیم دبی رودخانه بر 
سطح مقطع عرضی محاسبه می‌شود. بزرگی سرعت موثر نسبت به 
سرعت متوسط باعث شده گرادیان سرعت در رودخانه افزايش يافته و 
ضریب پراکندگی موّثر به همان نسبت بزرگ‌تر باشد. سینگ (,5108 
3 با ایجاد تغییراتی در ۸1۱8 مدل انتقال - انتشار اصلاح شده 
(۷۸۲۲۳۲) را برای رودخانه‌های دارای نواحی ماندابی توسعه داد. 
همچنین او با بررسی رودخانة میسوری» نشان داد نتایج مدل ۸۲ با 
جایگزین کردن سرعت و ضریب پراکندگی موثر در مدل به جای 
پارامترهای اصلی بهبود می‌اد 
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68 /20ع5]۵ -4 
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شکل ۱- طرح شماتیک از نواحی ماندابی در یک رودخانه (الف) نما از بالا (ب) نمای از کنار رودخانه 
۷۲ 56000۵ (0) 2۳00 ۲۱۵۲۷ هام (ه) فعصمم )صفصمهای عم مدز ماما -1 عساعز۲ 


سینگ (2008 ,510820) سه مدل انتقال آلاینده در رودخانه‌ها 
شامل مدل نگهداشت موقت (153) مدل انتقال - انتشار اصلاح 
شده (0۷۸۵) و مدل انتقال- انتشار کسری (۳8۸) را برای 
تعدادی از رودخانه‌های ایالات متحده آمریکا به کار برد و نتیحه 
گرفت که معادلة انتقال- انتشار اصلاح شده در مقایسه با دو مدل 
دیگر دارای پارامترهای کمتر و به لحاظ کاربردی ساده تر است. 
همچنین او نشان داد این مدل در مدل‌سازی رودخانه‌های دارای 
نواحی ماندابی بسیار تواناتر می‌باشد. 

طرح مسئله به عنوان یک مدل یک بعدی باعث شده بسیاری از 
عوامل دخیل در انتقال مواد محلول در رودخانه ساده شوند و یا در نظر 
گرفته نشوند. در شرایطی که داده‌های مشاهداتی برای یک رودخانه 
محدود باشد و يا مدل‌های پيجيدةٌ دوبعدی و سه‌بعدی نمی‌توانند به 
طور قابل اعتماد مورد صحت‌سنجی قرار گيرند» از سرعت و ضریب 
پراکندگی موّثر استفاده می‌شود. عموما در مدل‌های یک بعدی نواحی 
نگهداشت تعریف نمی‌شوند و اکثر مدل‌های یک بعدی بر اساس 
رابطه‌ی پایه‌ی انتقال - انتشار توسعه یافته‌اند. ارزش عملی این روش 
در این واقعیت نهفته است که اثر مناطق ماندابی در استفاده از 
معادلات یک بعدی سنت- ونانت و انتقال- انتشار اعمال شده و در 
مدل‌های توسعه یافته شده یک بعدی بکار گرفته شود. برای اینکه اثر 
نواحی ماندابی در معادلات دیده شود. لازم است مقدار 1 برای 
رودخانة مورد نظر تعیین گردد. اين پارامتر در نهایت از طریق سرعت 
(یا دبی) و ضریب پراکندگی موّثر وارد مدل می‌شود. رودخانة اروند 
علاوه بر مرزی بودن تنها ورودی آب شیرین به خلیج‌فارس محسوب 
می‌شود همچنین از آنجایی که اراضی و نخیلات اطراف رودخانه 


وابسته به کیفیت آب رودخانه هستند بنابراین کیفیت آب این رودخانه 
از اهمیت بالایی برخوردار است. با این حال به علت اینکه اطلاعات 
سطح مقطع عرضی و وضعیت بستر این رودخانه به طور دقیق در 
دسترس نیست و یکی از علل استفاده از مدل‌های یک بعدی همین 
موضوع می‌باشد» برای سنجش صحت روش پیشنهادی» از این 
رودخانه بهره گرفته شد. پارامترهای فاکتور پراکندگی (2) و توان 
پراکندگی (0) همراه با 17 از طریق واسنجی مدل 1۷۲116۳11 برای 
یک ایستگاه در اروندرود با الگوریتم بهینه‌سازی کلونی مورچگان 
تعیین می‌شود و برای دو ایستگاه دیگر مورد صحت‌سنجی قرار 
می‌گیرد. 


مواد و روش‌ها 
منطقةه موردمطالعه 

رودخانه اروند با طول تقریبی ۱۹۳ کیلومتر یکی از بزرگ‌ترین 
رودخانه‌های دجله و فرات در القرنه" نزدیک شهر بصره در عراق به 
هم می‌رسند و در ۳۳ کیلومتر شهر بصره رودخانه کارون به آن 
کیلومتری در مرز کشور ایران و عراق به خلیج‌فارس می‌ریزد (شکل 
۲ عرض این رودخانه از ۲۵۰ متر در پایین‌دست دجله و فرات تا 
بیش از ۷۰۰ متر در مصب رودخانه می‌رسد عرض این رودخانه در 
ورودی به خلیج‌فارس ۱۲۵۰ متر نیز بیان شده است. 
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ارو ندرود 
حاوش اد اعد 


(۳۱۰۰ 


۱۰ 


استگاه های اندازه گیری ‌ 
عصه‌تاد؟ 


رودخانه عم 
۷۰41 


شکل ۲ محد وده مطالعاتی و موقعیت رودخانه اروند 
۸۷۵۵/۵ 0 م۱0۵0 20 مع2 ت56 -2 ۲121۲6۵ 


ارتفاع کف این رودخانه در حدود ۰,۷متر در طول ۶۲ کیلومتر تا 
منطقه بصره کاهش می‌یابد. عمق رودخانه بین ۸,۵ متر تا ۲۴ متر در 
فصول پرآبی و کم‌آبی بیان شده است :2015 ,.۵1 6 طداآنتده) 
(2017 ,۸1-۸620 قسمت پایین حوضه دجله - فرات نشان دهنده 
بر کین مه خالای قر خایرماته و شرت آنببا اس که کت 
شرایط طبیعی ۲۰ هزار کیلومتر مربع را دربر می‌گیرد ,۱(۳۲۲) 
(2001. اصلی ترین منبع تأمين آب شیرین در آن منطقه رودخانه 
اروند است و پایدار ماندن بسیاری از نخیلات وابسته به کیفیت آب 
این رودخانه است. مهم ترین منبع کاهش کیفیت آب رودخانه نفوذ 
آب شور از سمت خلیج فارس به رودخانه است که بررسی عوامل اثر 
گذار بر برآورد کیفیت آب در اثر نفوذ آب شور را اجتناب ناپذیر کرده 
است (2023 .41 ۶ ۷]۵۵2261). 


مدل‌سازی 


در این پژوهش باتوجه‌به اطلاعات از وضعیت شوری و این که 


1 - ۱۷۷ ۵11-0 


(2017 ,.01 6۶ شبیه‌سازی جریان و کیفیت آب و انتقال آلودگی (در 
این مطالعه شوری) رودخانه توسط مدل یک بعدی انجام شد. 
مدل‌های زیادی برای تحلیل و شبیه‌سازی کیفیت و انتقال آلودگی در 
رودخانه‌ها توسعه داده شده است. 1۷]11۳11 یک ابزار یک بعدی با 
هیدرودینامیک قوی برای شبیه سازی جریان و کیفیت مصب‌ها 
روشخانفها و کال ها ام گه توسظ مقس فیدر ولیک نارق 
ایجاد و توسعه داده شده است (2016 ,111). به منظور برآورد غلظت 
شوری در رودخانه مورد مطالعه» دو ماژول هیدرودینامیک (11) و 
انتقال- انتشار (۸۲) مورد استفاده قرار گرفته است. قدم اول در 
مدل‌سازی عددی هر پدیده‌ای شناخت روابط حاکم بر آن پدیده است. 
ارم سیل دروبیک فتامل غمی آبا ورسرجت وریان از 
ورودی‌های مدل کیفیت برای پیش‌بینی شوری در مصب ارون درود 
می‌باشد. معادلات حاکم بر جریان رودخانه معادلات یک بعدی سنت 
ونانت است (2018 ,.۵1 6 صناتع]۷). 


(1۲۳۱۲) مانطتافما متانیجیل ۲ طوتمونز[ -2 
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در اين روابط ۸ سطح مقطع جریان ]4 4 دبی جریان ورودی 
جانبی در واحد طول [" ]2]» ۵ دبی جریان [ 1۳1 ه شتاب 
گرانش زمین [*۷1۲ ٩:‏ شیب خط انرژی, 50 شیب بستر " عمق 
آب [نلا» «مکان [آ] و ۶ زمان [7] می‌باشد. معادلة انتقال- انتشار 
برای پیش‌بینی و نحوةٌ انتقال مادهٌ آلاینده بسط داده شده که بر 
اساس معادله بقای جرم و قانون اول فیک" به‌دست‌آمده است. این 
معادله بیانگر دو مکانیزم انتغال جرم شامل انتقال ماده آلاینده توسط 
جریان و پراکندگی بر اساس اختلاف غلظت و سرعت است که در 
ماژول ۸ به‌صورت یک‌بعدی با فرض اختلاط کامل استفاده می‌شود 
(20058 ,20۲۵)). 
20060 5 (0)۸ 

0 0 

» غلظت آلاینده در شبکة ِ ۲ ۳ رودخانه ۸,173 
سطح مقطع جریان [آ]؛ ۵ دبی رودخانه []1]» 1 ضریب 
پراکندگی طولی [*1۲]» 6 ضریب زوال آلاینده (بعد آن تابعی از 
درجة واکنش می‌باشد. برای واکنش درجه صفر ["۷117*1]» برای 
واکنش درجه یک دارای بعد [17]» برای واکنش درجه‌دو "۲" ]1۷/]] 
[)» 5 ترم منبع [11771]» ۶ فاصله از محل تخلية آلاينده [1] و ؛ 
زمان [1] است. در مدل 1۷۲110۳11 حل عددی روابط دیفرانسیلی (۲) 
تا (۴) با روش تفاضل محدود به شکل ضمنی صورت می‌گیرد. 
انتخاب گام زمانی و مکانی تا حد زیادی برای ارضا کردن شرایط 
پایداری معادلة (۴) اهمیت دارد. شرط پایداری مدل 1۷1۲1۳11 در 
عبارت (۵) قابل مفاهده ات به علت اینکه جریان آب غیرمانتکار و 
غیریکنواخت در نظر گرفته شده است» شرط پایداری در هر زمان و هر 
مکان باید برقرار باشد. 


+45 3 ۸۲6 ۸5 )۴( 


98۸ 1 ۵ 


۸99 

قاتن اه از گام مات زار کم مان [۳۳ و تسش 
متوسط جریان در رودخانه [111] است. به علت اينکه ضریب 
پراکندگی میزان زیادی به گرادیان مکانی سرعت وابسته است. در 
مدل 1161311 این ضریب به صورت زیر تعریف می‌شود: 

0 "رام -< 

در این رابطه 7 سرعت متوسط جریان رودخانه [1/1] 2 فاکتور 
پراکندگی [127] و ط توان پراکندگی می‌باشند. در حالتی که ط برابر 


1- ۳10169 ۷ 


صفر باشد ضریب پراکندگی در طول ثابت و برابر فاکتور پراکندگی 
می‌شود و درصورتی که 0 عددی غیر صفر باشد» ضریب پراکندگی در 
محاسبات مدل وابسته به میزان سرعت بوده و در طول رودخانه و در 
زمان به علت غیریکنواخت بودن جربان متفیر خواهد بود. 

به منظور مقايسة و بررسی بهبود نتایج در مدل 1۷۲116۳511 ابتدا 
مدل بدون تغییرات پیشنهادی در رودخانة اروند اجرا شد. با توجه به 
اطلاعات موجود برای سال ۲۰۱۴ (شکل ۲ ابتدا سری زمانی دبی و 
غلظت مواد محلول ورودی به رودخانه به عنوان شرایط مرزی 
بالادست وارد مدل شد (2017 ,۵000 280 ۷]0۳260). با توجه 
به جداول جزرومد در دهانه خروجی اروند» تراز سطح آب دریا به 
عنوان شرایط مرزی پایین دست رودخانه در نظر گرفته شد ,5ع1201) 
(2014 و همراه با غلظت متوسط شوری خلیج فارس به مدل معرفی 
شد. قابل ذکر است شرایط مرزی پایین دست به صورتی تعریف شد 
که شوری خلیج فارس امکان بازگشت به شبکه محاسباتی در رودخانه 
را داشته باشد. همچنین» رودخانه کارون نیز به عنوان ورودی به 
شبکه‌ی رودخانه به صورت نقطه‌ای با جریان و غلظت ثابت تعریف 
شد (2022 ,.۵1 ۵۶ تععهاد۷۲0). 

پارامترهای هیدرودینامیک و کیفیت رودخانه پارامترهای مهم و 
اثرگذار در پیش‌بینی شرایط جریان و غلظت در رودخانه هستند. در 
این تحقیق ضریب زبری مانینگ ۰/۰۱۷ تا ۰/۰۳۲ و ضریب 
پراکندگی ۳۸۵ مترری بر ها ید مطالسانت و رارق ان 
1 201 .]60 61 وت ثِ«۳ :2019 .۵ 6 1261" 
,.1 61 ۱۷0۵/۵2۵1 2022). که در این پژوهش نیز از همین ضرایب 
برای شبیه سازی حالت پایه مدل ۱1116211 استفاده شد. 

به علاوه از اثر باد صرف نظر شد و ضریب واکنش با توجه به این 
که شوری آلاینده‌ای زوال‌ناپذیر است در شبیه‌سازی صفر (0< ) در 
نظر گرفته شد (2015 ,.۵1 2 1>2002). با عنایت به رابطه (۶) در 
بدا فاکتور پراکندگی بگونه ای تعریف شد که در طول رودخانه بدون 
تغییرباشد (0< ). بنایراین پس از مدل‌سازی هیدرودینامیک جریان 
(۳11) نتایج آن به مدل کیفیت (۸1) وارد شد. در نتیجه سری 
زمانی غلظت در طول رودخانه قابل استخراج است. 

حال با توجه به هدف مطالعه برای بهبود نتایج» تخمین بهتر 
تغییرات غلظت و جبران کمبود داده در مدل‌سازی باید تغییراتی در 
پارامترهای ورودی اعمال شود. تأثیر نواحی ماندابی در رودخانه از 
طریق تعریف پارامتر 7 و اثر آن بر سری زمانی جریان ورودی به 
رودخانه و ضریب پراکندگی دیده می‌شود. ضریب پراکندگی و سرعت 
موّثر به صورت عبارات (۷) و (۸) وارد مدل می‌شود. 
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شکل ۳- جریان ورودی به رودخانه اروند از بالادست در سال ۲۰۱ 
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در روابط ( ضریب پراکندگی []۲2]], ,1 ضریب پراکندگی 
موثر [" [*1] و 1 نسبت مساحت مقطخ عرضی ناحية ماندابی به 
مقطع عرضی رودخانه است. برای رودخانة اروند سرعت موثر از طریق 
تقسیم سری زمانی دبی ورودی به عبارت (1-7) تولید می‌شود. وارد 
کردن ضریب پراکندگی موثر به مدل ۱11211 کمی پیچیده‌تر به 
نظر می‌رسد. با توجه به وابستگی این ضریب به گرادیان مکانی 
سرعت. 7 از طریق فاکتور پراکندگی در مدل وارد می‌شود. پارامترهای 
و 7 در مدل 1۷۲116۳11 برای ایستگاه 56۳20 با فاصله ۱۲۳ 
کیلومتری از الفرنه توسط الگوریتم فرا ابتکاری کلونی مورچگان 
واسنجی شده سپس تغییرات اعمال شده در دو ایستگاه دیگر در طول 
رودخانه با توجه به داده‌های موجود اندازه گیری در سال ۲۰۱۴ صحت 
سنجی می‌شوند. ارتباط غیرمستقیم بین مورچه ها از طریق مسیرهای 
فرومونی آنها را قادر می سازد کوتاه ترین مسیرها را بین لانه و منابع 
غذایی خود بيابند. این ویژگی الهام بخش تعریف الگوریتمی جمعیت 
محور است که می توانند راه حل‌های تقریبی برای مسائل بهینه 
سازی پیچیده پیدا کنند. جزء اصلی الگوریتم های ۸0۵ مدل 
فرمونی است که برای نمونه برداری احتمالی از فضای جستجو 
استفاده می شود (2019 ,ع21ان۹ 40ههد مع1(02). تاکنون روش‌های 
مختلفی برای واسنجی مدل‌های انتقال آلاینده استفاده شده است 
(2003 .01 61 0600561 :1998 ,امعلصیت۳۴ :2019 ,.اه 6 ۴۵92 : 
7 .0 1 0۳۳2۵( با این حال محققین موفق به انتخاب بهترین 


روش واسنجی نشده‌اند. در این پژوهش از یکی از روش‌های فرا 
ابتکاری مبتنی بر هوش جمعی استفاده شده است که در براورد 
پارامترها در یک مجموعه یکتا توانمند است (2005 محت۱ظ : 
6 ,.01 61 1۱(0۲120). 

پیاده‌سازی الگوریتم کلونی مورچه در چند مرحله صورت می‌گیرد. 
گام اول ایجاد مورچه‌های اولیه. ارزیابی» مرتب‌سازی و قراردادن آن‌ها 
در بایگانی می‌باشد. در اين مقاله ۵ مورچه اولیه ایجاد شده که هریک 
۳ مقدار برای فاکتور پراکندگی (0) و توان پراکندگی (0) و نسبت 
مساحت مقطع عرضی ناحية ماندابی به رودخانه () در بر دارند. در 
گام دوم باید ضرایب وزنی :۱۷ و احتمالات :7 محاسبه شود. این دو 
ضریب را می‌توان با عبارات مختلفی محاسبه نمود که رایج‌ترین روش 
محاسبةٌ آن‌ها به صورت زیر می‌باشد: 
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در اين روابط / شماره مورچه اولیه موردنیاز» ‏ تعداد مورچه‌های 
اولیه» 4 ضریب کنترل تابع می‌باشد. در مرحله سوم مدل احتمالی 
مربوط به هرکدام از متغیرهای تصادفی ساخته می‌شود. برای ایجاد 
توزیع نرمال نیازمند میانگین و واریانس داریم که میانگین عدد 
تصادفی انتخاب شده در نظر گرفته می‌شود و واریانس از فرمول (۱۱) 
قابل محاسبه است. در نهایت برای استخراج مدل احتمالی توزیع‌های 
نرمال با ساختار فرمول (۱۲) ترکیب می‌شوند: 


منتظری و همکاران. بررسی اثر رفتار نواحی ماندابی در مدل‌های یک بعدی هیدرودینامیک و انتقال مواد محلول در رودخانه‌ها ۰ ۶۶۹ 
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در این عبارات 7 شمارنده متغیر تصادفی» 7 شمارنده مورچه ۸ 
مقدار متفیر تصادفی» یا ضریب کتترل تابع است. در نهایت با 
نمونه‌گیری تصادفی از هر متغیر و تولید مورچه‌های ثانویه. مورچه‌های 
جدید و موارد موجود در بایگانی ادغام شده و سپس آن‌ها را با توجه به 
انم هن میطه ی کب کروهه اسفای افنافق ختق می دسا 
بروزرسانی بهترین پاسخ یافته شده تا این لحظه می‌توان الگوریتم را 
در صورت نیاز تکرار نمود. در شکل ۴ فلوچارتی از نحوه‌ی مدل‌سازی 
مدل 1۷۲116۳11 در دوحالت پایه و اصلاح شده با هدف مقایسه سری 
زمانی غلظت خروجی از این دو مدل ارائه شده است. 


واسنجی مدل ۷۲116۳11[ اصلاح شده برای توابع هدف در 
ایستگاه ٩620‏ در فاصلة تقریباً ۰ کیلومتری از مصب اروندرود 


توننط لگوریتم کلونی مورچگان انجام ده است در این پژوهش 
برای واسنجی مدل. شاخص 18۳ میان داده‌های غلظت 


۱۲۱04:1160 1 
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تغییرات زماتی شوری 


مشاهداتی و غلظت‌های شبیه‌سازی شده توسط 1۷۲116۳11 به عنوان 
تابع هدف تعریف می‌شود. برای پیاده سازی الگوریتم کلونی مورچگان 
برای شروع ۵ مورچه اولیه انتخاب و ۱۰۰ تکرار در نظر گرفته شده 
است. روند پیشرفت در الگوریتم بهینه‌سازی در شکل ۵ قابل مشاهده 
است. مقادیر پارامترهای برآورد شده به ترتیب ۰/۱۶۸ < ( ۲۷۳/۴ < 
 <2 ۰۹۴‏ می‌باشد. 

افزایش پارامتر 1 باعث شده سرعت آب و ضریب پراکندگی در 
معادلات سنت- ونانت و انتقال- انتشار افزايش یابد. افزايش سرعت 
موجب کاهش زمان حرکت آب در رودخانه یا به اصطلاح 172۷61 
6 می‌شود. علی‌رغم افزايش ضریب پراکندگی به واسطة تأثیر 
ناحيةّ ماندابی» کوتاه شدن زمان حرکت آب در رودخانه باعث شده تا 
اثرگذاری این پارامتر بر سری زمانی غلظت شبیه‌سازی یا مشاهداتی 
سری زمانی غلظت بالاتر رود. تغییرات حاصل از اعمال اثر ناحية 
ماندابی بر نتایج شبیه سازی توسط معادلات یک بعدی سنت - ونانت 
و انتقال < انتشار در ایستگاه ٩6۳27‏ در شکل ۶ به خوبی قابل 
مشاهده است. در اين شکل داده‌های مشاهده شده با شبیه‌سازی 
متوسط غلظت ماهانه در اروندرود در دو سناریوی با و بدون اثر دادن 
نواحی ماندابی با یکدیگر مقایسه شده‌اند. 


۷111611 )61/۸51( 
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تغسرات زعاتی شوری 
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شکل 2- فلوچارت شبیه‌سازی و مقایسه نتایج دو مدل پایه و اصلاح ده 
۱۵۵066 ۵011160 مره متععول ۲۵ ۵۶ وااناعع۳ عمط ۵۶ دموزتومصی 0صه )سمل 10۱۲ متاهاسصنه -4 مباعز۲ 
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شکل ۵- تغییرات شاخص بهینه‌سازی با پیشرفت در الگوریتم کلونی مورچگان 
صتمعاه تجمای اصه از معمهمان از (ربا۵۲۲۵) ممتاه‌س عزدا0۵ ص مععصفل -5 ۲1۵0۲۵ 
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شکل 7- سری زمانی غلظت شوری برای ماه‌های مختلف برای ایستگاه 56027 
5۵۳۱ )۵ عط)جمصه احمعهللنل ۲۵۲ ممتاحجععصی زاتصاهه اه ممنهو م1 -6 عاع۲1 


همان‌طور که از شکل ۶ مشسخص است استفاده از ضریب 
پراکندگی و سرعت موّثر در معادلات. با افزایش ارتفاع بيشینه غلظت 
مشاهداتی و افزایش شیب در شاخه‌های پایین رونده و بالارونده باعث 
بهبود نتایج مدل شده است. این تغییرات به این معناست که شبیه 
سازی توسط مدل یک بعدی اصلاح شده به داده‌های مشاهداتی 
غلظت نزدیک‌تر شده است. بهبود تخمین غلظت در شاخة بالارونده و 


خطای جذر میانگین مربعات۲ (۲۳6۴) و شاخص نش ساتکلیف؟ 
(5۲) در جدول ۱ قابل مشاهده است. افزایش دقت نتایج مدل با 
م0611 001612110 -1 
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منتظری و همکاران. بررسی اثر رفتار نواحی ماندابی در مدل‌های یک بعدی هیدرودینامیک و انتقال مواد محلول در رودخانه‌ها ۰ ۶۷۱ 


افزوده شدن پارامتر 7 کاملاً مشهود است و به‌عنوان‌مشال خطای 
آماری 1۷1517 از 2/۲۱ 006 به ۰/۴۳ 001 کاهش می‌یابد همچنین 
شاخص 51 نیز به طور قابل ملاحظه‌ای به عدد یک نزدیک 
می‌شود. بنابراین نتایج شبیه‌سازی مدل با پارامترهای بهینه شده در 
این ایستگاه (86۳2۳) یعنی متغیر در نظر گرفتن فاکتور پراکندگی در 
طول رودخانه با توجه به سرعت جریان و اثر دادن نواحی ماندابی در 
باعث افزايش دقت مدل‌سازی شده است. 

در این بخش نتایج برای دو ایستگاه ۳۵۷ و ۲۷660 با فواصل 
۰ و ۱۵۰ کیلومتری از مبداً رودخانه مورد صحت سنجی قرار 
می‌گیرد. جانمایی این دو ایستگاه در شکل ۲ مشسخص شده است. 
نتایج سری زمانی غلظت شوری در ایستگاه‌های ۳۵۷۷ و 2/6۵0] به 
ترتیب در شکل ۷ و ۸ قابل مشاهده است. در اين دو ایستگاه پپس از 
اعمال ضرایب واسنجی شده در ایستگاه ٩6۳20‏ نتایج شبیه‌سازی 
مدل اصللاح شده به داده‌های مشاهداتی غلظت نزدیک‌تر شده است. 
با توجه به این که در این دو ایستگاه کلیات سری زمانی به خوبی 
مدل شده است با این حال بیشترین خطا متعلق به بیشینه غلظت در 
سری زمانی می‌باشد. این اختلاف ممکن است ناشی از خطا در ثبت 


داده‌های اندازه‌گیری و يا حتی ورود زه آب‌های کشاورزی از مزارع و 
نخیلات حاشیه رودخانه باشد که عموماً تخمین و اندازه‌گیری آن‌ها 
امری مشکل است. اما مدل بهبود یافته با دقت بهتری شرایط ماهانه 
غلظت را پیش‌بینی می‌کند که تخمین آن در شرایط متفاوت 
جریان‌های ورودی به رودخانه و حتی بالا آمدن سطح دریاها و تغییر 
در دامنه جزرومد برای استفاده کشاورزی و حتی شرب اهمیت قابل 
توجهی دارد. 

به طور کلی. غلظت شبیه‌سازی شده و اندازه گیری در ۱۲ ماه از 
سال در حالت اعمال ضرایب بهینه شده ضریب همبستگی () را در 
فاصله ۱۵۰ کیلومتر از بالادست از ۰/۸۶ به ۰/۹۷ می‌رساند و ريشه 
میانگین مربعات خطا (115۳) ۱/۲۷ 001 بهبود می‌یابد. در ایستگاه 
دیگر یعنی ۱۲۳ کیلومتری از بالادست اروندرود ضریب همبستگی از 
۳ به ۰/۹۸ رسیده است و با توجه به خطای 18۷۲5۳ مدل اصلاح 
هیا خسای کرش عطق ماهاکه را کیت ناهارمه 
سنجی دو ایستگاه ذکر شده با سه شاخص خطای آماری در جدول ۲ 


محاسبه شده است. 
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نوم نماد 


نتیجه گیری 


این پژوهش با هدف بهبود شبیه سازی انتقال واه محلول در 
رودخانه‌های تحت سلطة فرایند پرکندگی یک چارچوب برای اصلاح 
مدل‌های یک بعدی مرسوم به منظور افزایش دقت مدل‌سازی و 
نایک کرفیخ ان به شرایط واقعی, ارائثه داده است. این چارچوب 
فاکتوری اصلاحی مبتنی بر نسبت سطح مقطع عرضی نواحی ماندابی 
به سطح مقطع عرضی رودخانه را برای اصلاح سرعت جریان و 
وت وگن بر خی وک ایغ دا ارف > آلاان 
اعمال می‌کند. در اين مطالعه از مدل یک بعدی 1۷۲11۳11 برای 
رودخانه مرزی و جزرومدی اروند استفاده شد. در نهایت با به‌کارگیری 
تئوری فراابتکاری بهینه‌سازی الگوریتم مورچگان ضرایب مربوط به 
نواحی ماندابی و متغیرهای پراکندگی بهینه شدند و تغییرات به مدل 
معرفی شد. مقایسه نتایج شبیه‌سازی در دو حالت یعنی اعمال کردن 
ضرایب بهینه و حالت پایه با مقادیر مشاهداتی موجود. حاکی از دفعت 
بالای شبیه‌سازی کیفیت آب اروندرود در مدل اصلاح شده است. در 
مرحله بهینه‌سازی با مقایسه نتایج شبیه‌سازی با داده‌های مشاهداتی 
ضرایب ۰/۱۶۸ 2 1 ۲۷۲/۴ < ۰ ۰/۹۴ 2 تخمین دقیق‌تری از 


کیفیت آب برآورد کردند و در مدل اصلاح شده نسبت به مدل پایه 
خطای کمتری وجود داشت به‌عنوان‌منال خطای آماری 18151 پیش 
از ۵۰ درصد کاهش یافت و درنهایت نتایج شبیه‌سازی در دو ایستگاه 
و ۳2۷7 نیز با اعمال ضرایب به طور قابل ملاحظه‌ای بهبود 
داشت. 

نتایج بررسی نشان داد به طور متوسط در حدود ۱۶ درصد از 
سطوح مقاطع عرضی رودخانه اروند را نواحی ماندابی تشکیل می‌دهند 
که بیشتر این نواحی شامل تل ماسه‌های بستر و ناهمواری‌های موجود 
در بدنةّ رودخانه می‌باشد نکته‌ای که جلب توجه می‌کند کانال‌های 
آبیاری جزرومدی در جوانب رودخانه است که تعداد آن‌ها با توجه به 
مشاهدات میدانی و تصاوبر ماهواره‌ای نیز کم نیست. ورودی این 
کانال‌ها در مجاورت رودخانه منجر به تغییرات و به‌هم‌ریختن جریان 
آب در رودخانه می‌شود و عموماً سرعت آب در دهانه ورودی این 
کانال‌ها کمتر است. اين عامل به طور مستقیم ناحیه ماندابی نیست؛ 
اما اثرات کاملاً مشابه در طول رودخانه دارد و عموما اثر آن در 
مدل‌های یک‌بعدی دیده نمی‌شوند. با مقایسه سه ایستگاه اندازه‌گیری 
غلظت در طول رودخانه می‌توان گفت در ایستگاه‌های ۲6۵۲ و 
0 اثر نواحی ماندابی همچون پیچ و خم رودخانه به وضوح 


منتظری و همکاران. بررسی اثر رفتار نواحی ماندابی در مدل‌های یک بعدی هیدرودینامیک و انتقال مواد محلول در رودخانه‌ها 


مشخص است و در مدل اصلاح شده دقت شبیه‌سازی با تغییرات 
بیشتری بهبود یافته؛ اما در ایستگاه ۳2۷ به علت نزدیک بودن به 


اوه 


دخیل در مدل مورد بررسی قرار گیرد. از نتایج این پژوهش می‌توان در 
جهت دقت بالاتر در شبیه‌سازی‌های یک بعدی کیفیت آب و نزدیک 


کردن مدل‌ها به شرایط واقعی در رودخانه‌ها بهره برد. 


کمتر است» همچنین به علت نزدیکی به آب شور خلیچ‌فارس 
کانال‌های آپیاری کمتری احداث شده است. به همین علت اختلاف 
زیادی در مدل پایه و اصلاح شده رخ نداده است. رودخانه اروندرود 
یکی از رودخانه‌های پراهمیت کشور است که با توجه به شرایط مرزی 
آن اطلاعات اندازه گیری شده محدودی در دسترس و به عنوان یک 
رودخانه نمونه در این پژوهش مورد مطالعه قرار گرفته است. در نهایت 
پیشنهاد می‌شود برای دیگر رودخانه‌ها با داده‌های اندازه گیری 
کامل‌تر اثر نواحی ماندابی در مدل‌های یک بعدی با تغییر ضرایب 
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۶ 1060101112011107 .(2013) هشب باه عک رل 060051 رفظ ,۷۷2۳ ر.ظ مصطزلن۷6 و.؟ رت۷۷۵28۵606 بن) ,۱66۲[مکز 
.5290-2۰ :(49)9 66۵۲۵۱ ۳۳۵۹0۱۲۵۵5 ۲۷۵/۵ .صصتوعتاو صتعصامصط 2 مصمله تهج 9۵061 عمه۲ماو کصعمتعصهه 

ط ومامصصهتهن ممهتماو اصمتقصها ۵ فلوررلفمه زوجم 02022 ۷۲ رتمطه۷]22 عک ویک متحصفصظ تطاوح‌طاصهمطکز 
۵ 6۱2665 ۲۳۱۵۵۱۵۷ ,209۵ ها عصنصنهادهه ۲۳۷۵۲۵ طاً ا0۲ح8ص 2‏ حمتاتاهم ۵۶ . عصتام0مصظ . لهعتامصصعط)هظ 
۰ ۱05://01۰0۲8/10.22055/1186.2022.6ظ ۰ (امعتادهاه طوتاع مظ طازه صفتوعظ ط) .وععاممر 

صد للم اصممتمامم ۵ ز0«مصه ممتله‌تک0 106۸ ,(2021) .1 رعاموظ ع راو رال .۲.۷۷ ,960 بلط مطملا یگ ب۷0 15 
۱۱۱۲0۵۱۱۱۱۱۵۱۱۲۵ 0۴ 0۱۱۵ 17۱۵۳۱۱۵۵۵0۳۵۵ 15۰ هه ۷۵اه طوبامتطلههته طاذ مصتصندع۱ مصتطه‌محط عصلونا و۲۱۷۵ 
۰ :(6)3 7 601۱ ۴۱۵۱۱6 7۵0 

۲ ۸۱-۸۲۵0 )92 مطا ما ممتفتططا ماتصتلهه وم ممصمیااکصد تع۳ز یتک ۵۶ امممروومووه 21مزم ۱ .(2022) ,هه باه ] 
1-۰ :(2)6 ۲ ۵621۱66 ۲۷۷۵۲۵۲ ۸۵/۱۵۵ 60 صداحاحتظ ۵۶ ]وعطاتمه رادم 

۵ ۲8۵۲۷۱6 1۱۵۵ 0۴ 17۱/9۵۸6۵5 17:۵ .(2019) ۲۱۱۰۷۰ مططاظ عک و.تا.] معطهط! .1 بلططناه ر.ک مطعتاماصهک ر.ظ.1.ظ بل1۷]2 
۰ ۷۷۵۳۱۵ ۲۸۲۲۲۴ 38۳ مطا ۵۶ مهمن0عع۲-۳۳۵6 ,۰۵۲ ول رز مومس امه ۵ عوم و70۵۵ 

۷۷۵۲۵۲۳ 07 1۳۵5۴0۳ ۵ ۳۵۳۵۵۱۵۵۵۵۵ 2018(۰) ۱۱۷۷۰ بصفاامط0و اک 9 رهم۱۷۲۵ وسا.ل مطلاته]( 
۰ 6 .1110011718 

-۵1 926 معط صا معبتامنتاه مهقاطاصهووه صفظ صا فعصقطه تمممتاتوهم‌صهمن 2017(۰) .۵.۲ ,000ها عک ,.]1۷. ۸-1 ,عصهط۱۷]۵ 
5)5(۰ 6660665 ۵۵/۱20 ۵ 0۱۵1 ۸۵۱۵۲۱ 1۳۵0۰ رجع۷ ۳ حاحتظ 

021 ۸۳۷۵۱0 عط ما ممتمطصد اتمه که ام0مصظ آهم‌تاممصمطاه]۷ .(2022) .و بهتع۷]۵ ی ,۷ بتتعط۱۷۲222 بیظ بلتهمهاهه]۷( 
0 .221-248 :(48)2 وعزه 0۱۱۵۱۲۵۵۱۱۵ ۵و0 ,۲۷۵۲ معط مصبامج وصقا ۵۶ طتصتاهه ج0 ]ععه ما 2040 ۲۱۷۵۲ 
۰ 10.22059/168.2022.334989 /09://001.012ظ .موه طفتاع‌ص طز صولوزعظ 

-علع11 ۵۶ ممتاز0ا20 صهه‌اومها 01 ۲/6۵8 2023(۰) ,۷ رتطاعع۱۷]۵9200 ی ریق بلت۷۵1 م۷ بتتعطه۱۷۲22۵ ری بلله2ه]ه۷]۵( 
10 11:6 0۴ 96۱۵76۵ ,۷۵۲ ۲ع0۳0ظ طدتضصاه لفط ۵۶ وممتاموع۲ لذمو قح ماج ده صتفوظ ‏ تعرانک . ومهتطام‌بظ 
1۰ :56 ۱۱۷۱۲0۲۱۱۱۵۸۲ 

ماصمصتهاوناه 10 وععاعقعظ ممتنهام202 صمتویاصا للم ,(2019) ۳۰۵۰ ,عصوظ ع .۷ بتطعلط ۷2۵ .۲ ,۷۵ مصل .1۲.1 ,عصنطل 
اک ۱۱۱۱۵۵ .حصقهصام1 ۷ بدلنافمتمعط ]۱۷ من ۵۶ عوده ظ :عاهعه معصقطه متقصتاه بم4صن اصممصممام۷ع0 تناها هه 
68-۰ :23 ۲۲۵۱۱۵86۱۱۵۱۲ 

مصتتهتتاوی 20 ممتفتتاصا القو عصتاام۷۵0 ,(2022) ,۰1۷۲ رقهانا ی ر.۲ ,۷۲۵۳۵2۵00 ویکاما م۷1۵1 بن).۷۲ بهفتا90 بلط بقتله۳۵۲ 
00۱۵ .ا8دهه متاصدانه ۱۷۷ معط صم مومع تقصمتاتعصها مه صا دملتقصهعو فعصقطه ماقصته میا ففصصتاام 
 0)2( ۰‏ 1۱8۱66۳۱۱8۵ 0 ۹62۱۱6۵ ۲۲۵۲۱۱۱۵ 

6 ۶ ۷۵۲1۵008 110 از راصتهمعصه تمامصصه ۳2 2019(۰) ۲.۲۰ ,زماوع۲۱ 6 .۱۳.1 ما60 رانا بتلل[مع06ظ ,]9.۷ مفصف 
09۳۱۱۵ .مهن فامما۱۷ صتقطه ۷مته۷ عصتویا فحصحعه القصصه ‏ تم عمط ت0حقصه فتطامو . لقط10صممتل-عه۵ 
۰ 71۳01.2019«طر.[/۱05://601.0۲8/10۰1016ظ . بجوه۵] ۳0 

۳ ۱۱۵0۵46 ۱۲۵۲۱۵۵۵۲۴ 50/9۵ ۸ 0۵1156(۰) 10۲۵8۵ ۵۳00 سمل ۱ ۳۵/۵۵۲ 00۵-0۵۹۵۵ .(1998) ایک رلفلصنن۳۴ 
۰ 0۵01021021 ون۲۲] ,عمتتعاص عطع ۵۶ اممصص‌ندمعنا دت۲] ,(4018 مناوع1 ,98 ,۵1 ۱۷) ۳۷۵۲۶ ۵10 5۳۵۵/15 

۲ 00119۵۱۵۲۷۱ ۵ 0۴ (۱-1۱6) 6۵۵0 ۵6/۵۷۵۲۵ظ 2022(۰) .1 متصحصصح ۷۲۷ 6 ,.۷۲ ,۱۷۲222۳ ,.۵.۷۲ ,22021 
۱ 01 20۹( ۱ ۱۵24 20۳ 8 29 11 2 
۰ -۲۱09://01.012/10.22067/1661.2022.45 

عصطولصانظ اهصملووم۲۶۵1 نان .ععزهانه متسه ۱ ععنا هه هی (2005) 1۰۳۲۰۵۲۰ بوزنمه 52۷ 

٩0۵۲226 208 ۱۵061 10۳ ۲‏ ما 10۶ وعامصصدتهم ۵۶ ممتاهانمالیه 60وده-ممجصرم]۱۷ 2001(۰) .1.5 ,ق۵۵08 کی .1.۷۷ ,560 
453-۰ :(27)6 ۲ و6۱ع۱ ۱و ۵۳۱۲۵۵6 ول .صمتوه‌ووتل 

0 ۵9۳۱۱۵ .صمفممو0-ص0۷6000ع ‏ مصتنع1۳۷ طاً فمصم2. اصقصعهاو ۵۶ اصمحصاهع1:6.  )2003(.‏ کنو رطوصنه 
66470۰ (2003) 33-9429 ۸05://001:0۲8/10۰1061/)2866(07ظ 29)6(۰ ۲ ووععاع۲ 

۵9۱۱۵ ,۳2500۲۲ مصل۲۷6 صا 20865 اصحصوهاگ ۵۶ اصمطصصاهع۱ 10 قلعم معط ممتنهم‌هم (2008) کنو رطعصنه 
1/)2866(0733-9437)2008(134:6)853(۰ 10۰106 /6۵۹://001.0۲8 34)6(۰ ز وو ما ععم هه وه ممنمعز۳] 
عصتلسامص) 962 فصن طنامو هه صوعع0 صعتصا 3 آ۷۵ ,(۸1۲1) فعاطه) ما طلحهتطلجه ۱۳2۵3 .(2014) .]۸ روعاماه 1[ 
۰ ۵۵۵۲۱۱۱۵۱۱( ۲۵0رد۸۷۵ ۱6 1۵ ۵۲۵8۳۵۵/۵۲ ۰(وقاطاه 1 حصهعت 5 11021 


18. 


19۰ 


20 


21. 


22 


23. 


24 


25. 


26 


27. 


280 


29. 


30 


31۰ 


32. 


33. 


34 


5 


36. 


37. 


38 


39. 


40 


41۰ 


42. 


منتظری و همکاران. بررسی اثر رفتار نواحی ماندابی در مدل‌های یک بعدی هیدرودینامیک و انتقال مواد محلول در رودخانه‌ها 


۴ 0۳۱۵ ۱۷۵۲۹۰ صا امد صمتاالآه گه صمتاه 106208 ۲۵۲ 0مطامج 0عموجه-)صمجصرم]۱۷ 2015(۰) 2 رعصهی1 ی ر.ل ,1082" 
۰ :(10) ۲ 14 و۱66۳ ۱و۱ ۳۲۱۱۷۱۱۵۱۱۱۵۵ 

۲۷۷۵۲۱۱۵۵ اه ۵۴ ماع مومع صح ۵۶ موتصصع :ولصهاطوعه۷ صهتصصدا2۵0۵۵ع۱۷ م1 .(2001) .۳۱ رل 
۰ ۵۱۵ ۷۵۲۵۵ (0۷۶)) ۲و۵( ۵0۱۵ ق۷۵۵0 ۳۱۱۱۵0 957۵و 

بک بظ۲2 عک ویسا نک باق( ر.ظ مصصان۱۷۱6 ,رل ,۷0100۷۲۵ و.] رت22606 ۷۷ ویکل.ل.9 ب290 یت رتقتام(امک .رون ,۷۷۵2۲0 
ماه مامم۱۷ مصلودا وتمامصدنهم امه 0۲2۵0)و-)صع1فصه صاً پراصتهاتهمهصنا وقعوفه 6۵ ۵۵1 50109278 ۵ .(2017) بر 
4۰ 1103:0110 5661۱66۰ ۳۲۵/۱۷۷۵۵۲ 1201086۰ احط1و 

عم 2 مه آهءتعمامطامممه مه صمتعتطط تماجه 2 2006(۰) .۲ رما عک مرگ باملاط بل م92۷0 .میک مطط ۷۷ 
,137-9 :(22)1 6۶۵۵۲61 0۵81۵ 0۴ ۵۵۵ ۱6۲۲۱۱۵۲۷۰ حعطاتها ربع۷ ۲ ۲ماجع ناخ ]۴28 عط ۵۶ طاتامحظ 

-مصصی معتقا و نارس تقها فصقو ۵ ۵۶۴ حمتاقم؟ فطع عصتام۱۷]۵ .(2012) .ظ رعطتصصصما۲ عک رو رمقت ول م۷۷20 ب.م6 بلا۷ 
637۰ :299 «و660108 ۷۵۲۲۸۵ .مصنطن ,رتحتادظ عصه1زم‌صماصداه :۰تهتااوع 106-009182]60) ٩260,‏ 

0105 90۱۷و صرح سمل بع۲۱ مطا ممع0عه متطفممتهاع مطا عصتامعتادع۷ظ1 .(2022) .۱۱ رلط‌طاه)مک ی .۷۲ بطه 722 
89-۰ :(75)2 ۲ ۱۷/۵۱۵۵۵۱۱۵ ۲۷۷۵۸۵۲ -۱۵۵۲ع1 از 0 ۱۵۵۱ 0۱۱6 ۴۵0۵۵2 .صمتاحتاصععصهع 

معصههه معهتماو اممتمصهتا 0۶ اممصصووعو۸۵ 2016(۰) .0 محصممط2اتا نا ی و.ل ر16۳/۵000۷/۹1 ره مط۱۷]2110 بب۷ مقاله722721۳ 
۵ص .مه طمدامتطلهع:ها عصمله ومممصعته ممتاحجهعومع مهم فصوتصهده‌مصظ ‏ مملفتهمرفلل 20۷۵00 هه 
794-1۰ :538 ۳۱۲0/0 
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